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Abstract 


This  report  presents  an  overview  of  the  software  developed  under  the 
common  high  performance  computing  scalable  software  initiative  (CHSSI),  v 

computational  fluid  dynamics  (CFD)  6  project.  Under  the  project,  a 
proven  zonal  Navier-Stokes  solver  was  rewritten  for  scalable  parallel 
performance  on  both  shared  memory  and  distributed  memory  high 
performance  computers.  At  the  same  time,  a  graphical  user  interface 
(GUI)  was  developed  to  help  the  user  set  up  the  problem,  provide  real¬ 
time  visualization,  and  execute  the  solver.  The  GUI  is  not  just  an  input 
interface  but  provides  an  environment  for  the  systematic,  coherent 
execution  of  the  solver,  thus  making  it  a  more  useful,  quicker  and  easier 
application  tool  for  engineers.  Also  part  of  the  CHSSI  project  is  a 
demonstration  of  the  developed  software  on  complex  applications  of 
interest  to  the  Department  of  Defense.  Results  from  computations  of  10 
brilliant  anti-armor  (BAT)  submunitions  simultaneously  ejecting  from  a 
single  Army  tactical  missile  and  a  guided  multiple  launch  rocket  system 
missile  are  discussed.  Experimental  data  were  available  for  comparison 
with  the  BAT  computations.  The  CFD  computations  and  the 
experimental  data  show  good  agreement  and  serve  as  validation  for  the 
accuracy  of  the  solver. 
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COMPUTATIONAL  FLUID  DYNAMICS  (CFD)  COMPUTATIONS  WITH  ZONAL 
N A VIER- STOKES  FLOW  SOLVER  (ZNSFLOW)  COMMON  HIGH  PERFORMANCE 
COMPUTING  SCALABLE  SOFTWARE  INITIATIVE  (CHSSI)  SOFTWARE 


1.  INTRODUCTION 

The  thrust  of  the  work  described  here  is  to  further  develop  an  existing  computational  fluid 
dynamics  (CFD)  code  and  make  it  more  accessible  for  engineers.  The  code  was  developed  as  part 
of  the  common  high  performance  computing  scalable  software  initiative  (CHSSI)  and  is  now  called 
the  zonal  Navier-Stokes  flow  solver  (ZNSFLOW).  ZNSFLOW  is  actually  a  suite  of  codes  that 
basically  includes  a  zonal  Navier-Stokes  solver  and  graphical  user  interface  (GUI)  environment  for 
problem  setup,  interactive  visualization,  and  solver  execution.  The  primary  goals  of  the 
ZNSFLOW  CHSSI  project  are  to 

1 .  Develop  a  scalable  version  of  a  zonal  Navier-Stokes  solver; 

2.  Add  features  to  the  ZNSFLOW  software  which  allow  general  applicability  and  ease  of 
use;  and 

3.  Demonstrate  the  design  utility  of  the  scalable  ZNSFLOW  software  by  solving  current 
Department  of  Defense  (DoD)  priority  viscous  flow  problems. 

In  keeping  with  these  goals,  the  report  gives  a  broad  overview  of  the  ZNSFLOW  CHSSI  project, 
the  ZNSFLOW  solver  and  its  capabilities,  as  well  as  the  GUI  environment.  Some  results  from 
test  cases  are  presented  to  demonstrate  recent  applications  of  ZNSFLOW. 

2.  THE  ZNSFLOW  SOLVER 

The  ZNSFLOW  solver  originally  was  known  as  F3D,  a  fully  vectorized  Formula  Translator 
(FORTRAN)  77  code  used  on  Cray  vector  computers  such  as  the  C90.[l,2]  During  the  CHSSI 
program,  the  code  has  been  rewritten  to  provide  scalable  performance  on  a  number  of  computer 
architectures.  Other  added  enhancements  include  dynamic  memory  allocation  and  highly 
optimized  cache  management.  Aside  from  the  performance  aspects,  the  solver  has  been  provided 
with  a  number  of  enhancements  to  make  it  more  user  friendly  and  capable  of  performing  flow 
field  computations  for  complex  configurations  of  interest  to  DoD.  The  solver  was  written  to 
operate  with  and  without  a  GUI  environment.  A  large  portion  of  the  effort  spent  on  the 
ZNSFLOW  CHSSI  project  went  toward  increasing  ease  of  use  and  general  applicability  of  the 
ZNSFLOW  solver.  The  distributed  interactive  computing  environment  (DICE)  GUI  allows  some 
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of  the  ZNSFLOW  solver’s  more  complex  features  to  be  easily  employed.  For  example,  the 
solver  allows  for  one-to-one  overlaps  between  grid  zones  in  any  direction.  The  GUI  makes  the 
setup  for  this  generalized  data  exchange  intuitive  and  also  provides  some  simple  error-checking 
capabilities  to  catch  mistakes  in  creating  the  input  file.  Many  of  the  boundary  conditions  are 
generalized  and  can  be  applied  to  any  surface  or  line. 

In  addition,  the  solver  can  perform  computations  with  the  chimera  composite  grid 
discretization  technique.  [3-5]  By  using  the  chimera  technique,  one  can  greatly  simplify  the  grid 
topology  and  grid  generation  for  very  complex  systems.  One  of  the  drawbacks  in  using  the 
chimera  technique  has  been  the  increased  complexity  and  corresponding  confusion  in  applying  a 
turbulence  model.  A  chimera  model  can  be  composed  of  multiple  zones,  with  each  zone  possibly 
having  a  unique  grid  topology.  Most  turbulence  models  have  specific  directional,  orientation,  or 
distance-related  requirements  for  correct  application.  For  a  complex  chimera  model,  applying  a 
turbulence  model  can  be  a  very  complex  process.  This  problem  has  been  addressed  in 
ZNSFLOW  by  installing  a  pointwise  turbulence  model  [6]  that  is  not  orientation  specific.  This 
greatly  simplifies  the  setup  of  the  turbulence  model.  Wall  location  information  is  supplied  when 
the  wall  boundary  conditions  are  set  by  the  user  through  the  GUI.  Currently,  a  Baldwin-Lomax 
[7]  and  a  one-equation  pointwise  turbulence  model  have  been  validated.  A  two-equation 
pointwise  turbulence  model  is  undergoing  testing. 

The  ZNSFLOW  software  has  been  targeted  to  operate  on  both  shared  memory  and 
distributed  memory  architectures.  In  order  to  fulfill  the  CHSSI  program  requirement  that  the 
ZNSFLOW  solver  be  scalable  on  applicable  computers,  it  was  decided  to  create  two  versions  of 
the  ZNSFLOW  solver,  with  one  version  optimized  to  operate  on  shared  memory  architectures 
and  the  other  optimized  to  operate  on  distributed  memory  architectures.  The  shared  memory 
version  of  the  solver  employs  loop-level  parallelism  that  has  highly  optimized  cache  management. 
The  distributed  memory  version  is  not  as  fully  developed  as  the  shared  memory  version.  The 
distributed  memory  solver  currently  uses  the  shared  memory  (SHMEM)  library  for  the  Cray  T3E 
and  Origin  2000  computers  and  the  message-passing  interface  (MPI)  library  for  the  IBM®  SP® 
computer.  Either  current  library  can  be  chosen  with  conditional  compiler  switches.  Once  the 
ZNSFLOW  software  is  complete,  the  differences  between  the  multiple  versions  of  the  solver 
should  be  transparent  to  the  user.  Both  versions  of  the  ZNSFLOW  solver  apply  the  same 
unsteady  Reynolds  averaged  thin  layer  Navier-Stokes  equations  to  compute  flow  field  solutions 
with  no  changes  in  the  time-tested  solution  algorithm. 
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2.1  Governing  Equations  and  Solution  Technique 


The  complete  set  of  time-dependent,  Reynolds-averaged,  thin  layer,  Navier-Stokes 
equations  is  solved  numerically  to  obtain  a  solution  to  this  problem.  The  numerical  technique 
used  is  an  implicit,  finite  difference  scheme.  Steady  state  calculations  are  made  to  numerically 
compute  the  flow  field. 

2.1.1  Governing  Equations 

The  complete  set  of  three-dimensional  (3D),  time-dependent,  generalized  geometry, 
Reynolds-averaged,  thin  layer,  Navier-Stokes  equations  for  general  spatial  coordinates  T|,  and  £ 
can  be  written  as  follows  (Pulliam  &  Steger  1982): 

9rq +  d^F +  9^6  +  9^  =  Re'1  d;S,  (1) 

in  which 

^  =  £(x,  y,  z,  t)  -  longitudinal  coordinate; 

T|  =  ri(x,  y,  z,  t)  -  circumferential  coordinate; 

£  =  £(x,  y,  z,  t)  -  nearly  normal  coordinate; 
x  =  t  -  time 

In  Equation  (1),  q  contains  the  dependent  variables  (density,  three  velocity  components,  and  the 
energy),  and  F,G,andH  are  flux  vectors.  The  thin  layer  approximation  is  used  here,  and  the 
viscous  terms  involving  velocity  gradients  in  both  the  longitudinal  and  circumferential  directions 
are  neglected.  The  viscous  terms  are  retained  in  the  normal  direction,  and  are  collected  into  the 
vector  §.  In  the  wake  or  the  base  region,  similar  viscous  terms  are  also  added  in  the  stream- wise 
direction,  For  computation  of  turbulent  flows,  the  turbulent  contributions  are  supplied 
through  an  algebraic  eddy  viscosity  turbulence  model  developed  by  Baldwin  and  Lomax  [7]  or  a 
pointwise  turbulence  model  [6]. 

2.1.2  Numerical  Technique 

The  implicit,  approximately  factored  scheme  for  the  thin  layer,  Navier-Stokes  equations 
using  central  differencing  in  the  T|  and  C,  directions  and  an  upwind  scheme  in  £  is  written  in  the 
following  form  [2]: 
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/  +  ibh5\{A+)n  +  ibh8fn  -  ibh  Re-16cJ-IMnJ  -  ib  D,\  f 
I  +  ibhS{ (A-)n  +  ibhS1,Bn-ibD,\  „  A Qn  = 

,tAt{5|[(F+)”-#;]+  S{[(F-)"  -F;]+  «,[g”  ■ -G_]+  -  H.]- Re-%[S"  -S.]j 


(2) 


in  which  h  =  At  or  (At)/2.  The  free-stream  fluxes  are  subtracted  from  the  governing  equation  to 
reduce  the  possibility  of  error  from  the  free-stream  solution  corrupting  the  convergedjolution. 
Here,  5  is  typically  a  three-point,  second  order,  accurate  central  difference  operator;  d  is  a 
midpoint  operator  used  with  the  viscous  terms;  and  the  operators  8%  and  are  backward  and 
forward  three-point  difference  operators.  The  flux  F  has  been  eigensplit  and  the  matrices 
A,  B,  C,  and  M  result  from  local  linearization  of  the  fluxes  about  the  previous  time  level.  Here,  J 
denotes  the  Jacobian  of  the  coordinate  transformation.  Dissipation  operators  De  and  Dj  are  used 
in  the  central  space  differencing  directions.  The  smoothing  terms  used  in  the  present  study  are  of 
the  form 


and 

in  which 


DeL  =  (Al)J 


-1 


e2$p(B)fiS  +  e45 


Di\n  =  (At)J-'[e2Sp(B)p5  +  2.5e45p(B)S] 


I,  J> 

1 77  J '» 


and  p(B)  is  the  true  spectral  radius  of  B.  The  idea  here  is  that  the  fourth  difference  will  be  tuned 
down  near  shocks  (e.g.,  as  (3  gets  large,  the  weight  on  the  fourth  difference  drops  down  while  the 

second  difference  tunes  up). 


2.2  Chimera  Composite  Grid  Scheme 

The  chimera  overset  grid  technique  greatly  adds  to  the  number  of  applications  to  which  the 
ZNSFLOW  solver  can  be  applied.  The  chimera  overset  grid  technique,  which  is  ideally  suited  to 
multi-body  problems,  [8-10]  involves  generating  independent  grids  about  each  body  and  then 
oversetting  them  onto  a  base  grid  to  form  the  complete  model.  This  procedure  reduces  a  complex 
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multi-body  problem  into  a  number  of  simpler  subproblems.  An  advantage  of  the  overset  grid 
technique  is  that  it  allows  computational  grids  to  be  obtained  for  each  body  component 
separately  and  thus  makes  the  grid  generation  process  easier.  Because  each  component  grid  is 
generated  independently,  portions  of  one  grid  may  lie  within  a  solid  boundary  contained  within 
another  grid.  Such  points  lie  outside  the  computational  domain  and  are  excluded  from  the 
solution  process.  Equation  2  has  been  modified  for  chimera  overset  grids  by  the  introduction  of 
the  flag  ib  to  achieve  just  that.  This  ib  array  accommodates  the  possibility  of  having  arbitrary 
holes  in  the  grid.  The  ib  array  is  defined  so  that  ib  =  1  at  normal  grid  points  and  ib  =  0  at  hole 
points.  Thus,  when  ib  =  1,  Equation  2  becomes  the  standard  scheme,  but  when  ib  =  0,  the 

A  A  ,  I  A  A, 

algorithm  reduces  to  A Qn  =  0  or  Qn  —  Qn ,  leaving  Q  unchanged  at  hole  points.  The  set  of 
grid  points  that  forms  the  border  between  the  hole  points  and  the  normal  field  points  is  called 
inter-grid  boundary  points.  These  points  are  updated  by  interpolating  the  solution  from  the 
overset  grid  that  created  the  hole.  Values  of  the  ib  array  and  the  interpolation  coefficients  needed 

for  this  update  are  provided  by  a  separate  algorithm.  [3] 

Figure  1  shows  an  example  where  the  parent  missile  grid  is  a  major  grid  and  the  brilliant 
anti-armor  (BAT)  submunition  grid  is  a  minor  grid.  The  submunition  grid  is  completely 
overlapped  by  the  missile  grid,  and  thus  its  outer  boundary  can  obtain  information  by 
interpolation  from  the  missile  grid.  Similar  data  transfer  or  communication  is  needed  from  the 
submunition  grid  to  the  missile  grid.  However,  a  natural  outer  boundary  that  overlaps  the 
submunition  grid  does  not  exist  for  the  missile  grid.  The  overset  grid  technique  creates  an 
artificial  boundary  or  a  hole  boundary  within  the  missile  grid  that  provides  the  required  path  for 
information  transfer  from  the  submunition  grid  to  the  missile  grid.  The  resulting  hole  region  is 
excluded  from  the  flow  field  solution  in  the  missile  grid. 


Missile(Major)  Domain  Submunition(Minor)  Domain 


Artificial  Boundary  in  Parent  Domain 


Figure  1.  Inter-grid  Communications. 
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3.  THE  DISTRIBUTED  INTERACTIVE  COMPUTING  ENVIRONMENT 

As  stated  earlier,  ZNSFLOW  is  a  suite  of  codes.  Part  of  that  suite  is  DICE.[11]  At  present, 
only  the  shared  memory  solver  has  been  integrated  into  DICE.  DICE  provides  a  GUI  that  allows 
a  user  to  create  an  input  file  for  the  ZNSFLOW  solver.  DICE  can  also  be  used  to  execute  the 
ZNSFLOW  solver  once  the  input  file  has  been  created.  In  addition,  once  the  solver  is  executing, 
DICE  can  provide  real-time  visualization  of  the  flow  field  as  it  is  being  computed.  Even  if  the 
execution  of  the  solver  was  initiated  from  a  previous  day,  DICE  allows  the  user  to  connect  to  the 
application  on  a  remote  computer  and  visually  monitor  its  progress  on  a  local  workstation.  DICE 
provides  a  number  of  options  for  visualizing  data.  The  user  can  choose  surface  contours, 
isosurfaces,  x-y  plots,  or  spreadsheets  to  display  the  data. 

Figure  2  shows  some  of  the  GUI  windows  that  a  user  may  access.  A  flow  field  visualization 
window  is  visible.  The  user  may  interactively  rotate  or  translate  the  object  in  the  window  to  view 
the  flow  field  from  any  position.  Beneath  the  visualization  window  is  the  boundary  condition 
setup  window.  To  the  right  of  the  boundary  condition  window  is  a  data  directory  window.  This 
allows  the  user  to  drag  and  drop  specific  solver-generated  data  to  DICE  utilities  such  as  the 
isosurface  plotter.  Farther  to  the  right  is  the  solver  execution  window.  More  controls  for  excuting 
the  solver  on  multiple  platforms  are  available.  To  the  right  of  the  visualization  window  is  the  main 
interface  from  which  all  the  other  windows  are  initiated. 

It  is  important  to  note  that  DICE  is  not  only  a  GUI  but  an  environment  that  includes  a 
heterogeneous  distributed  memory  system  called  network  distributed  global  memory  (NDGM).[12] 
NDGM  uses  a  client-server  approach  that  allows  separate  distributed  applications  to  access  a  single 
contiguous  data  buffer  that  may  span  the  memory  of  several  computers.  Data  within  the  buffer  are 
stored  and  accessed  using  the  standardized  hierarchical  data  format  (HDF).  Data  from  the  solver 
can  be  placed  into  the  buffer  and  accessed  in  real  time  through  the  DICE  GUI  by  the  user.  DICE 
has  proved  to  be  an  exceptional  computational  environment  for  high  performance  computing 
software  and  is  currently  used  to  support  several  codes  developed  under  CHSSI  in  different 
computational  technical  areas  (CTA). 

4.  ZNSFLOW  DEMONSTRATION  CASES 

Demonstration  cases  were  chosen  to  show  the  capabilities  of  the  ZNSFLOW  software. 

Both  of  the  demonstration  problems  require  viscous  Navier-Stokes  CFD  modeling  for  accurate 
flow  field  solutions.  Two  demonstration  cases  were  run  on  a  Silicon  Graphics  Origin  2000 
computer.  The  first  of  the  two  cases  to  be  discussed  is  the  guided  multiple  launch  rocket  system 
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(MLRS)  missile.  The  guided  MLRS  computational  model  is  built  to  answer  questions  about  the 
use  of  canards  to  perform  controlled  maneuvers  for  a  missile  with  wraparound  tail  fins.  A  second 
demonstration  case  shows  the  capability  of  ZNSFLOW  to  model  complex  multi-body  systems. 
Computational  models  were  built  for  computing  the  flow  field  about  ten  brilliant  antitank 
submunitions  being  ejected  from  an  Army  tactical  missile  (ATACM).  The  complexity  and 
uniqueness  of  this  type  of  multi-body  problem  result  from  the  aerodynamic  interference  of  the 
individual  components,  which  include  three-dimensional  (3-D)  shock-shock  interactions,  shock¬ 
boundary  layer  interactions,  and  highly  viscous-dominated  separated  flow  regions. 
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Figure  2.  DICE  GUI  Window. 


4.1  Computations  for  Guided  MLRS  Missile 

The  computations  will  hopefully  provide  insight  for  engineers  into  the  interaction  of 
canards  with  the  down-stream  wraparound  tail  fins.  Providing  control  for  a  missile  with 
wraparound  tail  fins  is  more  complex  than  with  normal  tail  fins.  The  curvature  of  the  wraparound 
fins  allows  for  easy  storage  by  folding  against  the  missile  body  while  in  the  launch  tube. 

Immediately  after  launch,  the  fins  unfold  to  stabilize  to  the  missile.  The  cylindrical  shape  of  the 
wraparound  fin  is  advantageous  for  packaging,  but  it  can  also  compromise  the  dynamic  stability  of 
the  missile.  Wraparound  fins  have  a  number  of  unique  aerodynamic  traits,  the  most  infamous  of 
which  is  the  roll  moment  that  they  generate;  this  may  change  in  sign  and  magnitude  during  the 
course  of  a  trajectory.  The  roll  moment  contributes  to  the  missile  spin  rate.  During  the  course  of 
flight  of  a  wraparound  fin  missile,  it  is  possible  for  the  missile’s  spin  rate  to  increase  or  decrease 
more  than  once.  In  addition,  the  direction  of  spin  may  change.  This  type  of  behavior  can  produce 
poor  flight  dynamics.  CFD  can  be  a  useful  tool  for  predicting  the  aerodynamics  of  wraparound  fin 
missiles.  [13,14]  The  information  gained  from  the  computations  will  hopefully  aid  in  a  successful 
design  of  the  guided  MLRS  and  future  missiles  equipped  with  wraparound  fins. 

Initial  computations  have  provided  interesting  information  about  the  guided  MLRS  missile 
flow  field.  Wind  tunnel  data  are  available  for  comparison.  However,  because  of  an  error  attributable 
to  miscommunication,  the  nose  and  canard  geometry  of  the  computational  model  vaiy  slightly  fiom 
that  of  the  wind  tunnel  model.  The  results  still  provide  insight  and  will  demonstiate  the  capability 
of  ZNSFLOW  for  providing  flow  field  solutions  for  this  configuration.  The  computations  have 
been  run  at  0°  angle  of  attack  at  velocities  of  Mach  1.6  and  Mach  2.2  and  at  10°  angle  of  attack  at 
Mach  1.6.  For  all  computations,  each  canard  has  a  deflection  of  10°.  A  large  computational  model, 
exceeding  24  million  grid  points,  was  made  for  flow  field  computations  of  the  missile  at  angle  of 
attack.  The  computations  demonstrated  the  ability  of  ZNSFLOW  software  to  handle  laige  data 
sets.  The  computational  models  used  for  the  0°  angle  of  attack  computations  exploited  symmetry 
and  were  one-fourth  the  size  of  the  computational  model  used  for  the  angle  of  attack  case. 

Figure  3  shows  a  ZNSFLOW-computed  solution  of  the  guided  MLRS  missile  at  Mach  1.6  at 
0°  angle  of  attack.  Figure  4  shows  a  ZNSFLOW-computed  solution  of  the  guided  MLRS  missile  at 
Mach  1.6  and  10°  angle  of  attack.  The  flow  field  changes  substantially  with  the  increased  angle  of 
attack.  Figures  3a  and  4a  show  pressure  contours  on  a  plane  1.4  calibers  from  the  nose.  This  plane 
is  just  aft  of  the  canards.  The  location  of  the  vortices  generated  by  the  canard  tips  can  be  seen  as 
small,  circular  low-pressure  regions  near  the  canard  tips.  Figures  3b  and  4b  are  3.7  calibers  fiom 
the  nose.  The  flow  field  at  0°  angle  of  attack  is  symmetrical,  but  the  flow  field  for  the  10°  angle  of 
attack  is  asymmetrical.  Most  noticeable  is  a  large  low-pressure  region  on  the  visible  side  of  the 
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body.  Since  the  missile  is  flying  at  10°  angle  of  attack,  the  deflected  canard  beneath  the  body 
directs  more  air  to  the  visible  side  of  the  body.  This  low-pressure  region  extends  to  the  rear  of  the 
missile  and  is  visible  in  Figure  4c,  which  is  approximately  14  calibers  from  the  nose  and  is  just  in 
front  of  the  tail  fins.  Pressure  contours  in  Figure  3c  are  at  the  same  location  as  in  Figure  4c.  Figure 
3c  again  shows  the  symmetrical  flow  field  at  0°  angle  of  attack  in  contrast  to  the  asymmetrical  flow 
field  generated  at  10°  angle  of  attack  shown  in  Figure  4c.  In  Figure  3c,  the  light  contour  shade 
between  the  dark  contours  near  the  body  indicates  the  locations  of  the  tail  fins.  The  dark  pressure 
contours  in  Figure  3c  indicate  that  the  position  of  the  canard  tip  vortices  is  actually  between  the 
wraparound  tail  fins. 


a.  1 .4  calibers  from  nosetip.  b.  3.7  calibers  from  nosetip.  c.  14.4  calibers  from  nosetip. 
Figure  3.  Normalized  Pressure  Contours  at  Mach  1 .6  and  0°  Angle  of  Attack. 


a.  1 .4  calibers  from  nosetip.  b.  3.7  calibers  from  nosetip.  c.  14.4  calibers  from  nosetip. 
Figure  4.  Normalized  Pressure  Contours  at  Mach  1.6  and  10°  Angle  of  Attack. 


Visualization  using  particle  traces  has  also  provided  some  insight  to  the  guided  MLRS  flow 
field.  Figure  5  shows  particle  traces  released  from  the  wakes  of  the  deflected  canards.  The 
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particle  traces  for  Figure  5a  were  generated  from  a  Mach  1 .6  flow  field  solution,  while  the  particle 
traces  for  Figure  5b  were  generated  from  a  Mach  2.2  solution.  Figure  5a  and  5b  indicate  that  the 
flow  fields  are  similar  at  Mach  1.6  and  Mach  2.2  at  0°  angle  of  attack.  An  interesting  note  is  that 
particles  released  at  the  base  of  the  canards  nearly  hit  the  base  of  the  tail  fins.  However,  particles 
released  at  the  canard  tips  are  caught  in  a  vortex  that  passes  between  the  fins.  The  particle  traces 
in  Figure  5c  are  an  indication  of  the  differences  in  the  flow  field  for  a  guided  MLRS  missile  at  0° 
and  1 0°  angle  of  attack.  The  particle  traces  for  Figure  5c  were  generated  from  a  Mach  1 .6  flow 
field  solution  at  10°  angle  of  attack.  The  traces  show  that  particles  released  from  the  canard 
wakes  are  swept  to  the  lee  side  or  upper  side  of  the  missile  body.  As  mentioned  earlier,  the 
canard  beneath  the  missile  deflects  more  air  flow  to  one  side  of  the  body,  creating  a  large 
difference  in  the  flow  fields  on  the  sides  of  the  body.  For  the  10°  angle  of  attack  case,  only  the 
canard  tip  on  the  far  side  of  the  body  generates  a  strong  vortex.  An  indication  of  this  vortex  is 
the  small  dark  circle  on  the  left  side  of  the  body  that  is  visible  in  Figure  4b. 

4.2  Computations  for  RAT  Submunitions  Ejecting  From  ATACM 

The  ATACM-BAT  multi-body  problem  involves  the  radial  dispensing  of  several  BAT 
submunitions  (see  Figure  6)  at  a  low  supersonic  speed.  This  case  is  ideally  suited  for  the  chimera 
overset  grid  technique  described  earlier.  The  chimera  scheme  allows  each  BAT  to  be  modeled 
with  its  own  simple  orthogonal  grid  as  seen  in  Figure  7.  The  trajectory  of  the  3-D  radial 
dispensing  submunitions  depends  on  the  initial  ejection  velocity.  The  flow  field  is  complex  and 
involves  3-D  shock-boundary  layer  interactions  and  ATACM-to-BAT  as  well  as  BAT-to-BAT 
interactions.  Detailed  experimental  or  theoretical  data  were  not  available  to  help  evaluate  the 
submunition  dispensing  phenomenon  for  the  entire  BAT  system,  and  thus  the  numerical  solution 
of  this  problem  was  initiated.  [15-17]  The  chimera  solution  procedure  was  successfully  used  to 
help  determine  the  aerodynamic  interference  effects.[16] 

For  a  set  of  wind  tunnel  experiments,  the  position  of  the  submunitions  was  set  in  order  to 
evaluate  flow  field  correction  factors  for  nonsymmetrical  dispensation  at  a  distance  near  and  far 
from  the  bay.  The  flow  field  correction  factors  are  used  in  six-degree-of-ffeedom  simulations  of 
BAT  dispensation  for  different  conditions.  CFD  computations  were  made  for  two  configurations: 
Configuration  A,  which  places  the  submunitions  relatively  close  to  the  missile  bay,  and  Configura¬ 
tion  B,  which  places  them  farther  away  from  the  turbulence  generated  by  the  missile  bay.  For  both 
Configurations  A  and  B,  there  is  equidistant  circumferential  spacing  for  each  submumtion  except 
one,  which  has  a  5°  offset.  The  submunition  with  the  circumferential  offset  is  located  at  approxi- 
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mately  the  1 1  o’clock  position.  Figure  8  provides  a  visual  reference  for  the  submunition  positions 
for  Configurations  A  and  B. 


a.  Particle  traces  for  Mach  1.6, 0°  angle  of  attack. 


b.  Particle  traces  for  Mach  2.2,  0°  angle  of  attack. 


c.  Particle  traces  for  Mach  1.6,  10°  angle  of  attack. 

Figure  5.  Particle  Traces  at  Various  Mach  Numbers  and  Angles  of  Attack. 


Surface  pressure  contours  are  shown  for  Configuration  A  in  Figure  9  and  for  Configuration 
B  in  Figure  10.  The  surface  pressures  on  the  Configuration  A  submunitions  reveal  much  stronger 
pressure  gradients  than  the  submunitions  in  Configuration  B.  Also,  surface  pressure  contours 
within  the  AT  ACM  missile  bay  are  somewhat  different  between  Configurations  A  and  B.  The 
stronger  pressure  gradients  on  the  Configuration  A  submunitions,  which  are  much  closer  to  the 
ATACM  missile  bay,  are  indicative  of  the  higher  pitching  moments  generated,  which  tend  to 
push  the  nose  of  the  BAT  submunitions  radially  inward  toward  the  ATACM  missile  bay.  Since 
the  computations  include  multiple  BAT  submunitions,  BAT-to-BAT  interactions  are  included. 
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These  interactions  are  critical  and  have  a  strong  effect  on  the  aerodynamic  forces  and  moments. 
The  normal  force  and  pitching  moment  coefficients  vary  between  the  submunitions,  indicating  the 
asymmetrical  nature  of  the  interacting  flow  field. 


Figure  6.  Diagram  of  the  Multi-body  System. 
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Figure  10.  Normalized  Surface  Pressure  Contours  for  Configuration  B. 

Some  experimental  data  [18]  were  available  for  comparison  with  the  computational  results 
of  Configuration  A.  Figure  1 1  provides  a  visual  reference  for  location  of  the  BATs  that  were  the 
source  of  the  experimental  data.  A  BAT  at  approximately  the  5  o’clock  position  was  equipped 
to  record  pressure  data.  Pressure  data  were  collected  on  the  side  of  the  BAT  closest  to  the 
ATACM  and  the  side  facing  away  from  the  ATACM.  On  either  side  of  the  BAT,  pressure  data 
were  taken  at  five  positions.  Unfortunately,  the  pressure  data  obtained  from  the  experiment  on 
the  side  of  the  BAT  facing  the  ATACM  do  not  appear  to  be  accurate.  However,  the  pressure 
coefficient  data  computed  from  the  CFD  solution  on  the  side  of  the  BAT  facing  the  ATACM  are 
plotted  in  Figure  12.  Figure  13  shows  a  comparison  between  the  pressure  coefficient  obtained 
from  experimental  and  CFD-calculated  data  on  the  side  of  the  BAT  facing  away  from  the 
ATACM.  Both  Figures  12  and  13  show  the  pressure  coefficient  as  a  function  of  the  length  of 
the  BAT  body  where  X/L  =  0  corresponds  to  the  BAT  nose  and  X/L  -  1  corresponds  to  the  end 
of  the  BAT  body.  Figure  13  shows  that  the  pressure  coefficient  computed  from  the  CFD 
solution  is  in  very  good  agreement  with  experimental  data.  The  CFD-computed  data  plotted  in 
Figures  12  and  13  provide  an  interesting  comparison  that  demonstrates  the  asymmetry  of  the 
flow  field  about  the  BAT  and  the  strong  influence  of  the  ATACM  proximity  to  the  BAT. 
Although  the  comparison  between  the  experimentally  obtained  and  CFD-computed  pressure 
coefficient  is  quite  good,  the  comparisons  between  experimentally  obtained  and  CFD-computed 
force  and  moments  indicate  that  some  flow  field  characteristics  may  not  be  captured  accurately 
by  the  CFD  solution. 
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Configuration  A 


BAT  located  17.60  fsi 
from  ATACM  center 
where  force  and  moment 
data  were  collected .  ^ 

■t 


Inner  BAT  surface  where 
pressure  coefficient  data 
were  collected. 


BAT  located  14.85  fsi  from 
ATACM  center  where  force  and 
moment  data  were  collected. 


Outer  BAT  surface  where 
pressure  coefficient  data 
were  collected . 


Figure  11.  Locations  Where  Experimental  Data  Were  Collected. 


Force  and  moment  data  were  collected  from  the  BATs  located  at  the  12  o  clock  and  6 
o’clock  positions.  The  BAT  at  the  12  o’clock  position  has  a  radial  distance  from  the  ATACM 
center  of  17.60  full  scale  inches  (fsi).  The  BAT  at  the  6  o’clock  position  has  a  radial  distance 
from  the  ATACM  center  of  14.85  fsi.  Figure  14  shows  both  the  experimental  data  and  the  data 
computed  from  the  CFD  flow  field  solution.  The  data  in  Figure  14  indicate  that  the  CFD- 
computed  data  match  the  experimental  data  of  the  BAT  17.60  fsi  from  the  ATACM  center  more 
closely  than  the  experimental  and  computed  data  of  the  BAT  14.85  fsi  from  the  ATACM  center. 
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The  data  for  the  normal  force,  CN,  are  in  good  agreement  for  the  BAT  17.60  fsi  from  the 
ATACM  center.  The  side  force,  CY,  data  appear  to  be  the  same  for  the  CFD-computed  side 
force  and  the  experimental  side  force.  This  is  somewhat  misleading  because  the  magnitude  of  the 
side  force  is  much  smaller  than  the  normal  force  and  pitching  moment,  Cmz  (coefficient  of 
moment  about  the  Z  axis).  The  relatively  small  side  force  is  a  good  indication  that  the  BATs  are 
not  likely  to  move  closer  together  when  being  ejected  from  the  ATACM  bay  at  0°  angle  of 
attack.  The  difference  -between  the  pitching  moment  for  experimental  data  and  CFD-computed 
data  is  less  for  the  BAT  farthest  from  the  ATACM.  This  seems  to  indicate  increased  difficulty 
in  computing  the  flow  field  for  the  ATACM-BAT  multi-body  problem  accurately  when  the 
BATs  are  almost  in  the  ATACM  bay. 


Figure  13.  Pressure  Coefficient  Versus  BAT  Length  for  BAT  Surface  Facing  Away  From  ATACM. 
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Figure  14.  Force  and  Moment  Coefficients  for  Configuration  A. 
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The  drag  coefficient  computed  from  the  CFD  solutions  compares  very  well  with  the 
measured  drag  coefficient.  Figure  15  shows  a  plot  for  the  drag  coefficient  of  the  same  BATs  that 
were  instrumented  to  obtain  the  force  and  moment  data  displayed  in  Figure  14.  In  the 
experiment,  each  BAT  was  mounted  on  a  “sting.”  The  stings  were  not  modeled  in  the  CFD 
computation.  The  total  drag  coefficient  was  obtained  from  a  force  measurement  of  the  BAT  with 
sting.  The  experimental  value  of  the  BAT  forebody  drag  was  estimated  by  taking  a  pressure 
measurement  near  the  BAT  base  and  using  it  to  estimate  the  base  drag  component  of  the  total 
drag  coefficient.  The  base  drag  component  was  then  subtracted  from  the  total  drag  to  obtain  the 
forebody  drag.  An  interesting  note  is  the  increase  in  drag  with  the  increased  distance  of  the  BAT 
from  the  AT  ACM  center. 
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Figure  15.  Drag  Coefficients  for  Configuration  A. 

5,  CONCLUDING  REMARKS 

A  broad  overview  of  the  software  developed  under  the  CHSSI  CFD-6  project  has  been 
presented.  The  scalable  Navier-Stokes  solver  executed  through  the  interactive  computing 
environment,  DICE,  provides  enginneers  with  a  fast  and  comprehensive  CFD  computation  and 
analysis  tool  for  complex  configurations  that  require  large  computational  resources.  However,  it 
can  perform  computations  for  simple  cases  just  as  well.  The  software  allows  the  user  to 
perform,  monitor,  and  visualize  the  computations  on  large  high  performance  computers  without 
copying  the  computational  mesh  and  solution  to  their  local  workstation.  The  comprehensive 
interface  provides  control  for  every  aspect  of  the  computation.  It  was  also  demonstrated  that  the 
ZNSFLOW  software  provides  accurate  and  visually  informative  results  for  large  complex 
configurations  such  as  the  guided  MLRS  missile  and  BAT  dispersal  from  ATACM.  The 
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predictive  numerical  capability  documented  allows  for  accurate  computation  of  flow  fields  that 
capture  complex  aerodynamic  phenomena,  such  as  interference  effects,  required  for  the  improved 
design  and  modification  of  current  and  future  DoD  projects. 
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